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Abstract — Code-division multiple-access (CDMA) has the po- 
tential to support traffic sources with a wide range of quality 
of service (QoS) requirements. The traffic carrying capacity of 
CDMA channels under QoS constraints (such as delay guarantee) 
is, however, less well-understood. In this work, we propose a 
method based on stochastic network calculus and large system 
analysis to quantify the maximum traffic that can be carried 
by a multiuser CDMA network under the QoS constraints. At 
the physical layer, we have linear minimum-mean square error 
receivers and adaptive modulation and coding, while the channel 
service process is modeled by using a finite-state Markov chain. 
We study the impact of delay requirements, violation probability 
and the user load on the traffic carrying capacity under different 
signal strengths. A key insight provided by the numerical results 
is as to how much one has to back-off from capacity under the 
different delay requirements. 

I. Introduction 

The rapid growth of delay sensitive applications, such as 
VOIP and IP video, and the scarcity of the radio spectrum 
require the design of spectrally-efficient systems with qual- 
ity of service (QoS) support. Direct-sequence code-division 
multiple-access (DS-CDMA) |1| is an efficient technique that 
supports flexible scaling of the user population, bursty traffic 
sources, and a wide range of QoS requirements, such as delay 
guarantees. A fundamental question is: What is the traffic 
carrying capacity (throughput) of CDMA networks under the 
given delay constraints? Answering this question, however, is 
a challenging task because the existing physical layer channel 
models do not explicitly take into account the source burstiness 
and the QoS requirements. Furthermore, the time varying 
nature of the wireless channel may itself cause severe QoS 
violations. The above reasons call for a cross-layer approach 
that can model the link-layer user requirements (e.g. delay 
guarantee) and at the same time take into account the physical- 
layer techniques, such as error control coding |2| and multiuser 
detection 1 1 1 . Indeed, predicting the throughpuj^ regions for a 
given delay guarantee is seen as one of the grand challenges 
in multiuser networks |3|. 

Two front runner theories for cross-layer modeling of wire- 
less networks are the effective capacity [4], and the stochastic 
network calculus ||5), Q. The former has been used, e.g., for 
cross-layer modeling of single user multiple-input multiple- 
output (MIMO) channels with memory 17|. An attempt was 
also made in |j8) to use efficient bandwidth for cross-layer 

'Throughout the paper, the term delay refers to queuing delay and the 
throughput as the delay constrained throughput. 



modeling of wireless CDMA networks with linear minimum- 
mean square error (LMMSE) receivers. In [8], however, only 
memoryless channels were considered and the applicability 
of the proposed approach for quantifying the traffic carrying 
capacity of CDMA networks under delay constraints is not 
clear. Stochastic network calculus (NetCal) f9l, on the other 
hand, is a more general theory that has been recently applied 
for QoS analysis of wireless channels with memory Q, ||6j. 

In this paper, we propose a cross-layer approach to predict 
the traffic carrying capacity of multiuser CDMA networks 
under a given delay guarantee. We make use of moment gen- 
erating function (MGF) based stochastic NetCal that was pro- 
posed in fTO|, and further established and developed in 
In stochastic NetCal, a communication system is described 
by the arrivals at the ingress, the service it provides, and the 
departures at the egress. It enables easily the consideration of 
statistical independence between arrivals and system service. 
In contrast to the previous work, for example, in |8|, here 
the multiuser CDMA channel is considered to have memory. 
A significant amount of work has been done on link-layer 
channel modeling of such wireless channels (see p2) and 
the references therein), especially in the single-user setting. 
Most of this work builds on the finite-state Markov channel 
(FSMC) model, which has been successfully applied on the 
analysis of delay constrained single antenna |13| and MIMO 
I?] systems, to model the memory in the wireless channels. 
To the best of our knowledge, however, none of the works 
based on FSMC considers CDMA networks or uses the 
stochastic NetCal to analyze the delay constrained throughput 
of multiuser networks. 

The goals and major contributions of this work are as 
follows. We present a methodology to calculate the network 
delay constrained throughput of a multiuser DS-CDMA sys- 
tem that employs adaptive modulation and coding (AMC) 
at the transmitter, and LMMSE multiuser detection at the 
receiver The wireless channel is modeled via a FSMC and 
MGF based NetCal is used to compute the throughput of the 
system. The formulation is valid for any traffic arrival process 
for which the MGF or a bound on the MGF exists. In the 
analysis, we incorporate statistical independence of arriving 
traffic and the service provided by the channel. We show the 
impact of the delay guarantee, violation probability and the 
user load on the throughput under various conditions, such as 
signal strength and fading speed. 



II. The System 

The present paper considers a synchronous uplink CDMA 
channel with K users. Each mobile station (MS) is equipped 
with a single transmit antenna and the base station (BS) has a 
single receive antenna. Perfect channel state information (CSI) 
is assumed to be available at the BS, while the MSs have no 
direct access to CSI. At the physical layer, a discrete time 
signal model indexed by n — 1,2,... is considered. We let 
the number of active users Kn < K ?X any discrete time 
instant n be a random variable, modeling for example the 
bursty information sources. It is assumed that the BS knows 
the set of active users at all times. 

With the above assumptions, the discrete time received 
signal after matched filtering and sampling reads (T], | [T4| 

K 



1,2, 



(1) 



k=l 



where 6^, „ is the code symbol transmitted by the fcth user at 
time instant n. We let bk^n = if the user k is not active 
and fefe „ e C \ {0} otherwise. The code symbols of the active 
users are assumed to be independent with zero mean and unit 
variance. The length (the spreading factor) of the signature 
sequences 

{^k}k=i is M for all users. The random vector 
V ^ CN(0, ct^Ia/) represents samples of additive white Gaus- 
sian noise (AWGN), where Im denotes the identity matrix of 
size MxAI, and CN(/i, '4') stands for the circularly symmetric 
complex Gaussian (CSCG) distribution with mean fj, and 
covariance The channels between the MSs and the BS 
are assumed to be block fading (see, e.g., [15J ) with a fixed 
coherence time. Throughout the paper we assume that the code 
words transmitted by the users are shorter than the coherence 
time, and omit the block index for notational simplicity. The 
fading coefficients hk ^ CN(0,1) for all k — 1,...,K 
are considered to be independent and identically distributed 
(IID), but extension to unequal average received powers is 
straightforward. Random spreading |[T], | fT4| , |16| is assumed 
to be employed at the MSs, so that the signature sequences 
{sk}k=i are IID CSCG random vectors - CN(0, jjIm)- 
For future convenience, we define the user load of the system 
at time instant n as the ratio = Kn/M, and the average 
received signal to noise ratio (SNR) as SNRavg = 1/cr^. 

A. LMMSE Receiver for Multiuser CDMA 

Multiuser receivers mitigate, in addition to background 
noise, the interference between the users. It is well understood 
that the choice of multiuser detector has a significant effect on 
the system capacity (see, e.g., |[T|, p4| , 1 16|). The effect of 
delay constraints on the throughput of the CDMA network 
is, however, not yet well understood. Here we consider the 
impact of such effects on the throughput obtained by a mul- 
tiuser CDMA system equipped with the LMMSE receiver |lj. 
Extensions to other receivers are also possible, but they have 
been omitted from the present paper due to space constraints. 

Let us consider without loss of generality the detection of 
code symbol 6i.„ of the first user when the mobile stations k — 



1, . . . ,Kn < K are active. The LMMSE receiver forms the 
decision variable for the first user at time instant n as zi.„ = 
Ci„yt, where Ci „ e and the superscript {■)^ denotes 
for the conjugate transpose of the matrix. If the BS knows 
perfectly the channel coefficients {hk}^^i and the signature 
sequences {sk}^^i of all users, we can define 

Si,„ = [s2 ••• SKjeC*''"(^"^^\ (2) 
Di,„ = diag(|/i2p, . . . , \hKS) e c(^"-i»<(^"-i), (3) 

and the LMMSE detector for the first user becomes Q, |[14) 



(5) 



i + i/iiPstf'Mr^^.si 

The conditional output SINR for the LMMSE receiver reads 



ii,n{{sk}Li,{hk}Li) = |/iii's^Mr>i. 



(6) 



For the purpose of the cross-layer analysis, however, (|6]l is 
not a convenient starting point due to its dependence on the 
signature sequences and the channels of all active users. 

B. Large System Limit and the Decoupling Principle 

To simplify the physical layer model, we consider through- 
out the paper the large system limit in which both the number 
of active users Kn < K and the spreading factor M grow 
without bound with a fixed ratio a„ — Kn/M, for all 
n = 1, 2, .... In this asymptotic region, the multiuser CDMA 
channel ([T]l decouples under very general conditions to a set 
of single-user channels | [T4[ , p7| that do not depend on the 
signature sequences {sfe . For the case of LMMSE receiver 
considered in Section [Tl-A| the SINR after multiuser detection 
simplifies in the large system limit to p4) , p7) 

7l,n =Pl//3l,n, (7) 

where pi represents the instantaneous channel power in the 
decoupled single-user channel and is drawn according to the 
probability density function (PDF) 



fiPi) - 



Pi > 0. 



(8) 



The noise variance /3i_„ in the equivalent single-user system 
is the solution to the fixed point equation 



"dp. 



(9) 



/O P + Pl,n 

Since the interference becomes Gaussian in the large system 
limit pT) , the decision variables zi_„ are statistically equiv- 
alent to the received symbols of the single-user (su) channel 



n = l,2,. 



(10) 



where n™„ ^ CN(0, n — 1,2,..., are mutually 
independent samples of Gaussian noise. For later use, we 
denote „ = Epj^{7i^„} = l//?i,„ for the post-detection 
SINR (|7]l that is averaged over the channel gain pi . 

Remark 1: The fixed point equation ([9| is easy to solve 
iteratively and is always guaranteed to converge. Note also that 



/3i,„ has time dependence only through «„, so that in the large 
system limit only the size of the active user set matters — not 
which users are in it. For the purpose of cross-layer analysis, 
we shall therefore concentrate on the equivalent single-user 
channels defined by (|7]i - ( [TO| i, where {pk}^^i and {rfk ni^^i 
are all mutually independent random variables. 

C. Adaptive Coded Modulation for Multiuser CDMA 

Consider the single-user channel pO| , where the channel 
gain pi is fixed. Let the code word of the first user be 
bi — \t>i^i: . . . , hi_NiY G and denote the instantaneous 

received SNR over it 7^ = [71, „, . . . ,7i,ArJ^. Assuming the 
elements of bi belong to a discrete set A^, the maximum rate 
that can be reliably transmitted over ( [TO] i reads (see, e.g., |[T5)) 



jiNi) 



1 

(7i) = logzdXI) - log2(e) + ^ XI 9Mili,u), 



n=l 



as A^i — >■ cx). We denoted above 
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-|„ + y^(b_t)|2 



(11) 



dv, (12) 



and the integral is with respect to a complex variable v E C 
To introduce AMC to the system, we assume that: 

1) If the MS k is not active, it uses a code word = 0; 

2) For active users, the code word is selected from a 
sequence Ck.o,Ck,i, ■ ■ ■ ,Ck,L of random code book^ 
that have rates = i?o < -Ri < • • • < Rl < 00. 

All modulation sets Mi C C, ^ = 0,1,..., L, satisfy 
E |fofc.nP = 1, V6fe^„ € Ail, where E denotes expectation and 
Rq = "outage". For simplicity, the users k = 1,2,..., /f„ 
are assumed to be active at time instant n, and transmissions 
are initiated only at the beginnings of the fading blocks. With 
this assumption, the number of active users during a given 
fading block satisfies > K^i for all time instants n < n' . 

If the BS has knowledge of 7^, = [7^,1, . . . , ^k,NkV for all 
MSs, it can use ( fTT) to find the codes for the active user^ 
This, however, may consume too much uplink bandwidth and 
processing power at the BS. Here we we consider a simplified 
AMC scheme where, instead of using ([TTJ, the BS finds 



/* 



max 

!=0,1,...,A/ 



{l■■Rl<I']^S^k^)}, k^i, 



,Kn, (13) 



and feeds back the indexes to the active MSs before 

the data transmission starts. Each active MS then uses Ck,ii to 
transmit at rate over the CDMA channel with a vanishing 
probability of error as the code word lengths and coherence 
time grow large since (7fe,i) < lM'^\'lk) always holds. 

In theory, we would like to have a large set of code books 
with finely spaced rates to accurately adjust to the channel 

-For simplicity of notation, we assume tliat code words of any desired 
lengths can be picked from all code books Ck,i- 

^The users, e.g., inform their transmit payload before each fading block 
thi'ough a control channel, and the BS uses jvj (or l|6|) to compute 7^,. 



conditions. However, designing arbitrary rate codes is difficult 
and as the number of code books L grows, a reliable and very 
high rate feedback channel is needed for code book selection. 
Thus, the number of transmission modes L is in practice small. 

Remark 2: The authors in, e.g., [7], flE\ propose to sim- 
ulate specific AMC schemes in AWGN channel and use 
parameter fitting to derive "analytical" error probabilities for 
them. Such results are, however, heavily dependent on the 
chosen error control codes. We consider instead the modu- 
lation constrained capacity ( [TT] i that gives an upper bound for 
all coding schemes with the same modulation sets and code 
rates. With the new code designs 12J, we expect our results to 
provide a close description of modern wireless systems. 

III. Delay Constrained Throughput Analysis 

In this section we concentrate on the decoupled single-user 
channel ( [TOj l of the first user and omit both the user k and 
the time n indexes for notational convenience. The wireless 
channel service process in the block fading channel with AMC 
is modeled using an L-state finite state Markov chain (FSMC) 
where the block length (in seconds) is denoted by Tb. Each of 
the transmission modes is mapped to the corresponding FSMC 
state. Let Pc = [pi.v] be the L x L transition matrix of the 
FSMC. Assuming slow fading and a relatively small value of 
Tb, we have pi i' — for all |^ — > 1. The adjacent state 
transition probabilities are determined as |19| 

Njli + i) __ _ 

(14) 




; = i,2, 
;-i,2, 



.,L-l, 



where iV(7) denotes the level crossing rate (LCR) at the SNR 
value of 7. The remaining transition probabilities are given as 



Pi.i = 1 - P1.2 

PL,L = 1 -PL,L-1 
Pl,l = 1 -Pl,l-l -PIM 



1, / = 2,3,...,i-l. 



The level crossing rate reads for F > 1 19| 



N{V) 



/2^r^ / r 

7 V 7 



(15) 



(16) 



where /m — u/w is the maximum Doppler frequency of 
the channel, defined in terms of the vehicle speed v and 
the wavelength lu of the transmitted signal. The stationary 
probability tt; of the FSMC being in state / is given as 



r 



i+i 



7 



(17) 

The PDF of the instantaneous received SNR 7(7) of the 
single-user channel ( [TO] i is statistically equivalent to the large 
system post-detection SINR of ([T]), and given by 



1 / 7 
fil) = =exp -- 
7 V 7 



7>0, 



(18) 



where 7 = Epj{7i^„} — l//3i,„. We remind the reader that 
the average SNR 7 is in fact a function of time n = 1,2,..., 



c?x = inf < inf 

e>o T>o 



In 



,s-T)MsiO,s)^\ne <0 



(22) 



the noise variance at the receiver cr^, and the instantaneous 
user load through the fixed point equation (|9]). 

A. Main Result 

The moment generating function (MGF) of a stationary pro- 
cess X{t) is given by Mx{0,t) = E [e^-^^*)]. In the sequel, 
we denote Mx{0,t) = Mxi~d,t) for the parameters 9 > 0, 
t > 0. The discrete time arrivals and service of the channel are 
assumed to be independent stationary random processes given 
by the cumulative processes A{0,t) and S{s,t) respectively. 
For all real 9, the corresponding MGFs are M^(6',t) and 
Mg{9, t), respectively. Furthermore, let N\y denote the number 
of information bits per an upper layer data block and W the 
system bandwidth. Omitting again both the user and the time 
indexes for notational convenience, we can then write the 
number of data blocks transmitted in state I of the FSMC 
based channel service process as Ri = RiT^^W/N]^. 

Lemma 3.1: The MGF of the random process S{t) de- 
scribed by a homogeneous Markov chain with transition 
matrix and stationary state distribution vector tt, is given 
for 6* > and i > by 



TABLE I 
Transmission modes 



Ms{9,t)^7T{R{^9)P,Y- 
where 1 is column vector of ones and 
R(0) = diag(e^^S . . . , 



(19) 



OR 



(20) 

TTi] is given in ( fT?] ) 



The steady state vector tt = [tti tt2 
and Pc is obtained from ( [T4] l - ([TSj. 

Proof: We refer the interested reader to pO| . ■ 

We denote the traffic arrival rate by A while d| represents a 
bound with violation probability e on the delay. We can only 
provide a delay guarantee if d\ < 

Proposition 3.2: The delay constrained throughput Ad of a 
DS-CDMA system under a delay guarantee is given as 



Ad = a • max {A | < rf^}, 



where d\, assuming FIFO scheduling, is given in ( [22[ i at the 
top of the page and Ms{9,s) is given by ( [T9] l. 

Proof: We refer the interested reader to |f]T|. ■ 

Remark 3: The delay bound (|22]) is calculated using 
stochastic NetCal approach [111. MGF for a variety of arrival 
models is available in the literature flOl but finding the 
MGF of the service process Ms{9,t) is a challenging task. 
We address this challenge by making use of Lemma |3.1| to 
calculate Ms{6,t). Having obtained Ms{9,t), a stochastic 
bound on the delay (i| in (|22]) can be obtained using Chernoff 's 
bound, Boole's inequality and applying the technique in |TTj . 

We next use the result in this section to calculate the delay 
constrained throughput of a multiuser DS-CDMA system. 
Before closing this section, we present a definition which 
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will be used in the next section. The maximum carried traffic 
without the delay constraints in bps is given by 

L 

Cnn.^aWY,Rm- (23) 

1=1 

IV. Application and Numerical Results 

In the numerical results, we consider the AMC setup given 
in Table U The transmission modes follow the HIPERLAN/2 
specification (see, e.g., ||7|, |[T8|). The parameter 7/ in Table |l] 
denotes the SNR point in decibels (dBs) where the Itii mode is 
switched on. The physical layer of the system is parameterized 
by using a = 20 MHz channel. The base time unit is chosen 
to be Tb = 2 ms and a block transmission with fixed duration 
Tb is assumed. The size of the upper layer data blocks is 
fixed to iVb = 10 000 bits, but the number of information 
bits transmitted per fading block varies depending upon the 
selected AMC scheme. The rate matrix R(6') can be read from 
the Table |l] The infinite sum in the delay bound formula (|22]) 
is computed for the first 4000 units of time. 

The framework discussed earlier enables the derivation of 
delay constrained throughput of a various number of traffic 
sources whenever the MGF exists. In this work we use a 
periodic source for the discrete time model that generates 
arrival traffic. Such a traffic source with period t produces 
6 units of workload (data blocks) and its MGF reads |5| 



(21) MAi9,t) = e' 
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) i^'' 
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(24) 



where t > and 6 > 0. Here we set t = 1 and use the 
number of generated data blocks 6 to set the arrival rate A. 

Fig. [T]depicts throughput as a function of delay guarantee # 
for different delay bound violation probabilities e. One can 
observe the decrease in system throughput as the violation 
probability and / or the delay guarantee get(s) tighter 

The impact of the average received SNR (SNRavg = l/c^) 
on the throughput is studied in Fig.|2] As expected, increasing 
the average received SNR leads to an improvement in the 
throughput, as seen in the classical information theoretic 
results as well. For the considered system this is due to the fact 
that, the higher average received SNR increases the likelihood 
of a better channel (note that here this effect comes from (|7]i, 
(|9|, ( [T7] i and ( [T8] l) at any given time instant and, thus, allows 
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Fig. 2. /m = 50 Hz, a = K/M = 0.5, e = lO'^ 

for higher rate code books to be used at the transmitters on 
average. This can be clearly seen by observing the steady state 
vectors, e.g., for SNRavg = -2 dB and SNRavg = 4 dB 
which are [0.622,0.234,0.127,0.017,0.00044,1.43 x 10"^ 
and [0.25,0.184,0.261,0.203,0.096,0.01,3.991 x 10"^], re- 
spectively. 

Finally, Fig. [5] depicts the delay constrained capacity as 
a function of the maximum user load a — K/M. With 
the delay constraint, the system throughput behaves similarly 
to the information theoretic case (see, e.g., 1 ,16) Fig. 1]), 
but the optimum user load shifts towards zero when the 
delay guarantee gets tighter. The non-linear behavior of the 
throughput as a function of the maximum user load is due to 
the interplay between (|9]l, ( |2T| ) and ( |23| l. We also see that the 
higher the user load is, the further we have to back-off from 
the capacity to meet the delay guarantee. 

V. Conclusions and Future Work 

In this paper we formulate a method to find the delay 
constrained throughput of a multiuser DS-CDMA system with 
time varying channel. Channel memory is modeled by a finite 
state Markov chain. We present numerical results where we 
quantify the impact of increasing the delay guarantee and the 
user load on the delay constrained throughput for varying 
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signal strength. This work finds application in performance 
evaluation of wireless networks where the maximum through- 
put for a given delay guarantee and amount of resources is of 
interested. While in this work, we only used periodic source 
to illustrate the results, the same methodology can be applied 
to any traffic source with known MGF 
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